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Abstract: Infrared divergences in scattering amplitudes arise when a loop momentum I 
becomes collinear with a massless external momentum p. In gauge theories, it is known that 
the L-loop logarithm of a planar amplitude has much softer infrared singularities than the 
L-loop amplitude itself. We argue that planar amplitudes in N = 4 super- Yang-Mills theory 
enjoy softer than expected behavior as I || p already at the level of the integrand. Moreover, 
we conjecture that the four-point integrand can be uniquely determined, to any loop-order, by 
imposing the correct soft-behavior of the logarithm together with dual conformal invariance 
and dihedral symmetry. We use these simple criteria to determine explicit formulae for the 
four-point integrand through seven-loops, finding perfect agreement with previously known 
results through five-loops. As an input to this calculation, we enumerate all four-point dual 
conformally invariant (DCI) integrands through seven-loops, an analysis which is aided by 
several graph-theoretic theorems we prove about general DCI integrands at arbitrary loop- 
order. The six- and seven- loop amplitudes receive non-zero contributions from 229 and 1873 
individual DCI diagrams respectively. 
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1 Introduction 

It has long been appreciated that direct Feynman diagram calculations are often a very in- 
efficient way to study multi-loop scattering amplitudes, especially in a theory as simple as 
planar N = 4 super- Yang-Mills (SYM) theory. Rather, one usually aims first to obtain a 
representation of a desired amplitude as a linear combination of a (hopefully) small number 
of (hopefully) relatively simple integrals. To date there have been at least four essentially 
different technologies available for determining an integral representation for an amplitude. 
In this paper we introduce a new approach which is both conceptually simple and compu- 
tationally powerful, as we demonstrate by using it to determine the seven-loop four-particle 
integrand in SYM theory. 

1.1 Brief Review of Four Paths Towards Integral Representations 

The most conventional approach is to begin with a collection of integrands, make the ansatz 
that the amplitude should be a linear combination of those integrals, and then determine 
the coefficients of the integrands by matching various data — such as generalized unitarity 
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cuts [1-3] or leading singularities [4-8] — between the ansatz and the amplitude. Each 'data 
point' generates a linear equation on the coefficients, so for sufficiently many data points, and 
a sufficiently large set of linearly independent integrands, one can obtain a unique solution 
for the coefficients. Actually establishing the correctness of the resulting ansatz (i.e., ruling 
out the possibility of additional contributions to the amplitude which happen to vanish on 
all of the data points considered) usually requires careful analysis, for example via more 
complicated D- dimensional unitarity cuts. Details of this tried and true approach can be 
found in the review [9], and we have learned that its power has been exploited in as yet 
unpublished work to construct an integral representation for the four-particle amplitude at 
six-loops [10]. 

A very different approach has been introduced in [11, 12], in particular in [13], where 
a relation was presented which allows one (in principle) to write down the integrand of any 
desired multi-loop amplitude in SYM theory recursively in terms of lower-loop integrands 
(involving successively more particles). Note that we distinguish here between the integrand, 
which is a uniquely- defined rational function of internal and external kinematic data in any 
planar theory, and an integral representation, which is understood to be well-defined only 
modulo the addition of terms which integrate to zero. 

A third even more recent approach relies on the duality, or equivalence, between cer- 
tain correlation functions in SYM theory and maximally helicity violating (MHV) scattering 
amplitudes at the level of the integrand [14]. Aspects of this duality have been studied in 
several papers including [15-18]. For the special case of four particles the relevant correla- 
tion functions possess an extra hidden symmetry which has been exploited to provide a slick 
derivation of the four-particle integrand at three-loops in [19], and at higher loops in currently 
unpublished work [20]. 

A fourth approach makes use of the understanding of how soft- and collinear-infrared 
singularities factorize and exponentiate in planar gauge theory [21-34]. If we restrict for 
simplicity our attention to MHV amplitudes and denote by = A^ /Affl the ratio 

of the L-loop n-particle MHV (color-stripped, partial) amplitude to the corresponding tree 
amplitude, then exponentiation implies that in dimensional regularization with D = 4 — 2e 
we have 



even though the individual contributions on the left-hand side have stronger infrared (IR) 
singularities, Mn = 0(e~ 2L ). The general form of this statement is regulator-independent: 
for example, in the Higgs regulator of [35], the right-hand side of eqn. (1.1) would be C(log 2 m) 
while each individual M„ diverges as 0(log 2L m) as the mass is taken to zero — a fact 
which has been used to guide the construction of various integrands in the Higgs-regularized 
theory [36, 37]. 

The divergent behavior of equation (1.1) has often been used as an important consistency 
check on the correctness of integral representations [6, 7, 38-41]. Moreover it can be used to 
guide the construction of a representation for an L-loop amplitude by using the requirement 
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that all of the 0(e 2L ) through 0(e 3 ) poles must cancel as 'data' in the same sense as 
described above. For example, for the four-particle amplitude at three-loops, if one starts 

(3) 

with the one- and two-loop amplitudes as given, and makes the ansatz that should be a 
linear combination of the three-loop ladder and tennis court diagrams, then their coefficients 
can be uniquely fixed by requiring that eqn. (1.1) should hold at order A 3 . In this case, the 
resulting ansatz agrees with the correct three- loop amplitude found in [42] , though in general 
this approach is not guaranteed to work, since it cannot detect the contribution from any 
individual integral whose infrared behavior is too soft. 

1.2 Preview of a New Path Towards Integrands 

The last approach outlined above has the obvious disadvantage that it relies upon explicit 
evaluation of multi-loop integrals, a problem which — despite recent progress — remains ex- 
tremely difficult in general. In this paper we explore an approach which amounts essentially 
to imposing eqn. (1.1) at the level of the integrand, which allows us to work with simple 
rational functions instead of complicated polylogarithm functions. (Indeed it has been noted 
in [43] that soft limits can also be handled very easily at the level of the integrand, in par- 
ticular an (n + l)-particle integrand should reduce directly to an n-particle integrand in the 
soft limit p n —> 0. This statement is not generally true for integrals due to IR-divergences.) 

Infrared divergences arise when some loop momentum becomes collinear with an external 
momentum p a . In terms of the dual coordinates x a defined by p a = x a +i — x a this happens 
whenever some loop integration variable xa approaches the line connecting x a and x a +\ for 
some a (see [13, 44] for a detailed discussion). If we take a limit where both (x a — xa) 2 and 
(x a +i — xa) 2 are order e (but making sure that all other kinematic invariants take generic, 
non- vanishing values) then any integrand will have a 1/e 2 pole in the limit e — > 0. However, 
we claim that the integrand of the logarithm of any amplitude has only a 1/e pole at any 
loop-order L > 1. This phenomenon is clearly related to the fact of eqn. (1.1), but at the 
moment we present this claim only as an empirical observation which we have verified for the 
four-point integrand through seven loop-order. 

Assuming that this behavior holds for all L, then it obviously can be used as a powerful 
tool for constructing integrands of amplitudes. We simply need to postulate a suitable basis 
of integrands for any desired amplitude and try to find a linear combination for which the 
C(l/e 2 ) poles cancel. We call this procedure the soft-collinear bootstrap because it allows the 
determination of an L-loop integrand given all of the corresponding integrands at loop-order 
< L, which appear in the L-loop logarithm. (This approach has a similar flavor to the work 
of ref. [45], where it was argued that all tree- level amplitudes could be fully constrained by 
their collinear singularities, combined with the requirement of Yangian symmetry). 

The case of four particles is special because a very simple basis naturally suggests itself: 
the collection of integrands which are invariant under conformal transformations on the dual 
x variables, appropriately called dual conformally invariant (DCI) integrands [46, 47]. This 
class of integrands has already played an important role in guiding the construction of integral 
representations at four- [40] and five- loops [3]. We prove in Lemma 2 of appendix A that 
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every individual DCI integrand has a nonzero 1/e 2 singularity. We also have checked through 
seven- loops that the family of DCI integrands at each loop-order is linearly independent 
(if one imposes dihedral symmetry in the external particles, which is a symmetry of all 
superamplitudes; it is curious to note that imposing mere cyclic symmetry is sufficient to fix 
the amplitude through five- loops, but the full dihedral symmetry must be imposed in order 
to obtain the six- and seven- loop amplitudes). Moreover they remain linearly independent 
even after picking-off from each one just the residue of its 1/e 2 pole. Taken together, we 
are therefore led to conjecture that the four-point integrand can be uniquely determined, 
to any loop-order, by imposing the 0(1/ e) collinear behavior together with dual conformal 
invariance and dihedral symmetry. 

It is of course well-known that the four-point amplitude in SYM theory is trivial, in 
the sense that its logarithm is determined exactly by the anomalous dual conformal Ward 
identity [47] up to an overall multiplicative factor of the cusp anomalous dimension /(A) [48], 
whose value is in turn known exactly from other considerations [49]. However, since /(A) 
can obviously be computed by simply integrating the integrand, it is amusing to note that 
while (anomalous) dual conformal invariance of the logarithm of the amplitude is not powerful 
enough to fix /(A) alone, our results suggest that dual conformal invariance of the integrand 
is, when one also imposes the mild 0(1/ e) collinear behavior together with dihedral symmetry 
with respect to the external particles. Moreover despite the triviality of the amplitude, the 
integrand of its logarithm is of independent interest since it is evidently a nontrivial function 
which makes other remarkable appearances, for example as a Wilson loop expectation value 
in twistor space [50-53] and as (the square root of) a certain correlation function of stress 
tensor multiplets [14-18]. 

In section 2 we review important facts and conventions related to integrands in SYM 
theory and present our main conjectures (all of which we have verified through seven-loops). In 
section 3 we explain the classification of DCI diagrams, which we use as a basis for constructing 
the four-particle integrands. Our results are presented and discussed in section 4, and the 
reader may find more information at [54]. 

Throughout this paper we use the momentum twistor parameterization [55]. Conse- 
quently all momenta (both internal and external) are strictly four-dimensional and so all 
Gram determinant conditions are automatically satisfied. However we suspect that when 
expressed in terms of the x's, formulae for four-point integrands are actually valid for planar 
SYM theory in any number of dimensions, a property which is known to hold at least through 
four- loops [56]. 

All of the important ingredients in this paper, including the exponentiation of infrared 
divergences and dual conformal invariance, rely crucially on planarity, but it would clearly be 
of great interest if some variant of our approach might be applicable beyond the planar limit. 
Integral representations for the complete four-point amplitude in SYM theory, including all 
non-planar contributions, are known through four-loops (see for example [57]). 



- 4 - 



2 Conventions and Conjectures 



2.1 The Integrand 

The integrand of the four-point amplitude = /A^ in planar SYM theory is, at 

L-loop-order, a rational function involving only the scalar quantities 

x a b = (x a - Xb) 2 , with s = x 2 4 and t = X13 (2.1) 

formed from the dual coordinates xi,X2,xs,X4 associated with the external kinematics and 
the L loop integration variables xa, %b, For example at one and two- loops we have [58, 59] 

M[ 1 \ XA )= ™ 4 , (2.2) 

X1AX2AX3AX4A 

M?\x A ,x B ) = ^ + ^ + {A»B). (2.3) 

X±AXiBX 2 BX3AX3BX4AXAB X1BX2AX2BX3AX4AX4BXAB 

In order to avoid a little bit of clutter we omit throughout this paper the somewhat conven- 
tional overall prefactors of (—1/2) L and 1/L! — the latter from symmetrization with respect to 
the loop integration variables, and the former coming from a choice of how to normalize the 
measure of integration. In the interest of specificity, let us note that to compute amplitudes in 
the normalization convention of for example [59] our L-loop integrands should be integrated 
with the measure 

1 A L 1 f TT d 4 XA k 



n ■ (2-4) 



k=l 

The above examples exhibit three important general properties of integrands: 

1. full permutation symmetry in the L integration variables x A , xb, ■ ■ ■] 

2. dihedral symmetry D4 in the external variables xi, ■ ■ ■ , 24; 

3. and the absence of double poles. 

2.2 The Integrand of the Logarithm 

The above properties ensure that there is no ambiguity in defining the integrand of the 
logarithm of an amplitude. Taylor-expanding the left-hand side of eqn. (1.1) to seventh-order 
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in A, we have 

(logM„)( 2 )= M£)-\{MVf- 
(logM„)( 3 ) = - M^MW + ^MWf ; 

(logM„)( 4 ) = M W - mWmW - \{M^f + M( 2 )( M «) 2 - i(M«) 4 ; 
(logM„)( 5 ) = M( 5 ) - M^MW - m^mM + M( 3 )(M«) 2 + (M^M^ 
-M^{M^f + \{M^f- 

+ 2M^MP M W - M( 3 )(M«) 3 + \(M^f - \{M^f(M^f 

+ M( 2 )(M«) 4 -i(MW) 6 ; 

(logM n )( 7 ) = MP - mWmW - M^M® + M( 5 )(M«) 2 - A#>Af 

+ 2M( 4 )M( 2 )M« - M( 4 )(M«) 3 + (M( 3 )) 2 MW + M^\M^f 
-3M( 3 )M( 2 )(M«) 2 + M( 3 )(M«) 4 

- (M( 2 )) 3 M« + 2(M( 2 )) 2 (M«) 3 - M^{M^f + i(MW) 7 . 

For the two- and three- loop logarithms, these expressions are to be interpreted at the level of 
the integrand respectively as 



2 m-) 2 (2.5) 



(logM n ) (2) (^,x B ) = M^(x [A ,x B] ) - 1 -M£\x [a )m£\x b] 



M^{x A ,x B ) -^(mV){xa)mV>{xb) + mM(x b )mM{xa) 



(2.6) 



and 

(log M n ) W ( XAi XB , xc ) = {x [A , x B , x c\ ) - (x [[A , x B] )M« (x c] ) 

+ \m£Xx [a )M^(x b )m£)(x c] ) 
= MP{x A ,x B ,xc)-{MP{x A ,x B )MP{xc) (2.7) 
+ MP(x A ,xc)M£\ XB ) + M^(x b ,xc)mW(xa)' 



- ^M^ 1) (x A )M^ 1) (x i j)M^ 1) (xc) + 5 other permutations 

where we have used the notation [A ■ ■ ■ C] to denote the outer-symmetrization over indices 
A, . . . ,C (being careful to avoid over-counting when symmetrizing the product of separately 
symmetrized functions). Each of these expressions manifests the properties outlined at the 
end of section 2.1. 
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2.3 The Main Conjecture 

In order to phrase our main conjecture precisely, it is best to employ the momentum twistor 
variables of Hodges [55] . These parameterize the on-shell external momenta in terms of four- 
component momentum twistors Z l a , a = 1, . . . , n, i = 1, . . . , 4, which are related to the scalar 
momentum invariants according to 

x ab = (x b -x a 2 = y /. oc aa+1 66+1). 2.8 

(aa + l)(bb + 1) 

Here (as usual) all indices a, 6, . . . are understood modulo n and the bracket denotes the 
determinant 

(abed) = det (Z a Z b Z c Z d ) . (2.9) 

Finally the 'oc' in (2.8) indicates that we will not keep track of the two-brackets in the 
denominator. This is well-justified because all dependence on them drops out of any dual 
conformally invariant function of the x ab . 

Each off-shell loop integration variable xa is parameterized by the antisymmetric com- 
bination XA-i A xa 2 °f a P a h °f four-component momentum twistors. These can appear in 
brackets of the form 

x aA = (x a -x A ) 2 ^ (aa+lA 1 A 2 ), x A b = [xa - x B ) 2 -> (A\ A 2 B\ B 2 ), or etc., (2.10) 

where, as forewarned above, we omit the two-brackets completely since our intention is to 
only ever use these replacements in dual conformally formulae. 

We can probe the infrared structure of an amplitude at the level of the integrand by 
taking the limit as one of the loop integration variables approaches the line connecting x a 
and x a+ i for some a, as discussed in [13, 44]. Due to the symmetry of the integrand it is 
sufficient to consider the limit as xa approaches the line connecting x\ and x 2 . This can be 
accomplished by taking Za x to Z 2 and Za 2 to any generic point which lies on the hyperplane 
spanned by Z\,Z 2 ,Z^, i.e. 

Z Al ->Z 2 + 0(e), Z A2 ^aZ 1 +pZ 2 + 7 Z 3 +0(e). (2.11) 

In this limit, we clearly have 

x 1A = (12 A! A 2 ) oc e, x 2A = (23Ai A 2 ) oc e. (2.12) 

By 'generic' we mean that we don't want to accidentally choose Za 2 so that some other 
singularities are probed at the same time; in particular this means that we must take both a 
and 7 to be nonzero, for if a were zero then x%a = (34^4i A 2 ) would also vanish, while if 7 
were zero then x^a = (4 1 A\ A 2 ) would also vanish. 

It might be interesting to attempt to glean further information about integrands by 
probing these multi-collinear regions, but at the moment we content ourselves with simple 
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collinear limits. With this in mind, let us consider the generally safe choice of taking a,f5,j 
all to be 1: 

Z Al -> Z 2 + O(e), Z A2 ^ Z x + Z 2 + Z 3 +0{e). (2.13) 

In this region of the A\ integral, the statement that an integrand has at most an 0(1/ e)- 
divergence can be formalized as the conjecture: 

Main Conjecture. The integrand of the logarithm of the n-particle L-loop amplitude in 
planar SYM theory behaves as 0(1/ e) in the limit (2.13) for all L > 1. 

2.4 Two- and Three-Loop Examples 

The integrand of the 2-loop logarithm (log M^)^ is the sum of 5 terms, obtained by plugging 
eqns. (2.2) and (2.3) into eqn. (2.5). Only the three terms containing both x\a an d x 2 a m 
the denominator contribute to the 0(1/ e 2 ) pole in the limit (2.13). Combining these three 
terms over a common denominator yields 

n „W2)/ \ X13X24 (X1BX24X3A + X13X2BX4A- Xi 3 X 2 4X AB) , , x ,n 1 A\ 

(logM 4 y z >(xA,x B ) = hO l/e . (2.14) 

X1AX1BX2AX2BX3AX3BX4AX4BXAB 

Turning to momentum twistors, we can plug in 

324 = (2341) = -(1234), (2.15) 

x 3A = (3 4 A 1 A 2 ) = -(1 2 3 4) + 0(e) , (2.16) 

xi3 = (1234), (2.17) 

x 4A = -(1234) + 0(e), (2.18) 

so the nontrivial numerator factor becomes 

xibX24,x 3A + x 13 x 2 bX4A - X13X24X ab = (123 4) 2 (xib - x 2 b - x A b) + 0(e) , (2.19) 

but this in turn is easily seen to vanish as a consequence of 

x AB = (A x A 2 Bx B 2 ) = {21Bx B 2 ) + (23Bi B 2 ) = -x lB + x 2B ■ (2.20) 

The same cancellation of the 0(l/e 2 ) pole occurs at three- loops, but let us turn the 
statement around by pretending for a moment that we did not know the integrand, but were 
willing to make the ansatz that it should be a linear combination, 

mP(x a ,x b ,x c ) = a/{ 3) + b4 3) , (2.21) 

of the only two available dual conformally invariant integrands 

3 

4 3) = X13XU + sym _ ( 12 terms tota n ^ 2>22 j 

X\AX\BX\CX2CX3AX3BX3CX±BX abx AC 

if = X 2 13 X24X 2A + gym ^ (24 termg tQtal) _ (2 23) 

X1AXWX2BX2CX3AX3BX4AX ABX ACXBC 
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Assembling eqn. (2.7) leads to an expression for (\ogM 4 ) {3) (x A ,xb,xq) as a sum of 36 terms. 
If we isolate those terms containing both x\a and xia we find that the integrand vanishes at 
order 0(l/e 2 ) in the limit (2.13) only for a = b = 1, in accord with the known value for the 
three- loop integrand [42]. 

2.5 Integral Basis Conjectures for n = 4 

As exemplified in the previous paragraph, we can use our conjecture as a tool for determining 
an integrand, but only if we first identify a basis of integrands for constructing a suitable 
ansatz. For n > 4 there are two separate potential problems. First of all, because MHV 
amplitudes for n > 4 are maximally chiral, they cannot be represented by the parity-even 
quantities x a b alone; rather they must involve parity-sensitive generalizations of x a b such as, 

(A 1 A 2 2A) (2.24) 

which could perhaps be written u x+a" where were defined as one (of the two) complex 
points in space-time simultaneously light-like separated from all four momenta pi,...,p4. 
Secondly, the class of multi-loop integrands involving such factors tend to satisfy many 
integrand-level relations (that is, they generally are not linearly independent but rather form 
an over-complete basis of integrands). 

It is well-known that scalar box integrals provide a basis for integral representations 
of all one-loop amplitudes in SYM theory, modulo terms which integrate to zero in four 
dimensions, and it has recently been shown that certain chiral octagons provide a basis 
(though overcomplete) for all one- loop integrands [44] . The construction of a basis for planar 
two- loop integral representations has been carried out in [60], and it would be interesting to 
explore the implications of that work for the special case of the integrand of planar SYM 
theory. 

However for the special case of n = 4 we expect a very simple story: namely, that the 
integrand can be expressed as a linear combination of rational functions of the x a i,, each 
of which is invariant under dual conformal transformations [61] (in particular, under the 
simultaneous inversion x^ — > x^/x 2 of all external and internal dual variables). Indeed this 
property was used to guide the construction of the four- and five-loop amplitudes in [40] 
and [3], where such integrals were called 'pseudo-conformal'. 

Dual conformal invariance is broken by the dimensional regularization traditionally used 
to regulate the infrared divergences of loop- amplitudes, so one might expect that imposing 
dual conformal invariance at the level of the integrand is not exactly synonymous to imposing 
it at the level of the integrated amplitude. Indeed, it was observed in [62] that only a 
particular subclass of pseudo-conformal integrals actually appear, namely those which can be 
rendered infrared-finite by evaluating them in a particular scheme (off-shell regularization) 
which manifestly preserves dual conformal invariance. We use the abbreviation DCI to refer 
to an integrand which is both pseudo-conformal in the sense of [3] and finite when regulated 
off-shell as in [62]. 
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The observation of [62] through five-loops together with our analysis through seven-loops 
motivates us to make the following conjecture: 

Weak Basis Conjecture. At any loop-order, the integrand of the four-particle amplitude 
in planar SYM theory can be expressed as a unique linear combination of dual conformally 
invariant integrands involving only the quantities (x a — x^) 2 . 

Let us emphasize again that this conjecture makes two logically separate claims: first of all, 
that for any L, the integrand can be expressed in terms of this particular class of objects; 
and secondly, that this collection is linearly independent for any L. Also, we note that 
this conjecture is stronger than a corresponding statement about dual conformally invariant 
integrals. For example, if two quantities differed by something which integrates to zero, they 
would be linearly independent as integrands but not as integrals. 

We prove in Lemma 2 of appendix A that every DCI integrand has a non-vanishing 
C(l/e 2 ) singularity in the limit (2.13). Therefore, it is not possible for any integral to escape 
detection using our method of probing the collinear singularities of integrands. However, in 
order for the method to have maximal power, which is to say in order for it to be able to 
uniquely determine the coefficient of every single DCI integrand, it is necessary for us to make 
the slightly stronger conjecture: 

Strong Basis Conjecture. At any loop-order, the collection of DCI integrands remains 
linearly independent if we isolate from each one the residue of the 1/e 2 pole in the limit (2.13). 

We have verified that this is true through seven-loops by explicit calculation. 
3 A Basis for the Integrand 

In this section we outline the classification of the dual conformally invariant diagrams we use 
as a (conjectured) basis for constructing four-particle integrands. 

3.1 Pseudo-Conformal Diagrams 

Since all quantities constructed from the x a b's are automatically invariant under translations 
and rotations in x-space, the only nontrivial constraint from dual conformal invariance arises 
from imposing invariance under inversions, x« — > Xa/x 2 , which takes 

/ _ \2 . ( x a ~ x b)^ _ %ab ,4 d^X , * 

x ab — \ x a x b) ' 9 9 — ~~ 9 9 > a X r . „. , . Ky-'-) 

In order to check whether a given rational function of the x a bS is invariant under this transfor- 
mation, one simply has to check that all of the x 2 factors which accumulate as a consequence 
of eqn. (3.1) cancel out. This is tantamount to simply counting how many times a given index 
appears in the denominator and numerator. For example, for the first term of the two-loop 
integrand (2.3), 

^13^24 ,„ oA 



Figure 1. Dual conformally invariant diagrams through three- loops. We use the standard notation 
where a dotted line connecting two faces x a , Xb corresponds to a numerator factor x a b, but for clarity 
we omit from each diagram an overall factor of X13X24. These four diagrams correspond to the rational 
functions given in eqns. (2.2), (2.3), (2.22) and (2.23) respectively. 

we note that indices 1 and 3 each appear twice both upstairs and downstairs, while indices 
2 and 4 each appear once both upstairs and downstairs, so all of their associated x\ weights 
cancel-out. At the same time the indices A and B associated to the integration variables each 
appear four times downstairs, which is expected as their weights should cancel against the 
factors arising from transforming the measure d 4 XAd 4 XB- 

At this level of analysis we are looking simply at the naive transformation of integrands 
under (3.1), with no consideration given to the question of infrared divergences which arise 
when the integral is actually performed. Following [3] we therefore refer to any quantity 
invariant under (3.1) as pseudo-conformal. In figure 1 we review a standard diagrammatic 
notation for expressing the class of pseudo-conformal quantities of interest. (However at high 
loop-order we abandon this notation because the proliferation of numerator factors renders 
most diagrams illegible.) Two minor differences with respect to the standard notation are 
that: (1) for clarity we omit from each diagram an overall factor of X13X24 (we prove in 
Lemma 1 of appendix A that every four-point DCI diagram has such an overall factor); and 
(2) for each diagram we impose (as emphasized above) by construction that the resulting 
rational function of x a bS should have D4 x Sl symmetry and should be normalized so that 
the numerical factor in front of each distinct term in the sum is precisely 1. 

The translation between a diagram and its associated rational function of x a ft' s proceeds 
via the following simple steps. We label the four external faces xi, X2, x$, X4 in cyclic order 
and then label the L internal faces of an L-loop diagram as xa, x b, ■ ■ ■■ For example, let us 
choose to label the fourth diagram in figure 1 as: 



Each propagator bounded by faces x a and Xf, is associated as usual to a factor l/x ab , while 
each dotted line connecting faces x a and Xfj is associated to a numerator factor x a fy. (Note 




that, as in the third diagram in figure 1, there may be k > 1 dotted lines connecting the same 
two faces, in which case the numerator factor x a b is raised to the k-ih power.) This diagram 
has 10 propagators and 2 numerator factors, so according to the rules it gives the factor 



where we have restored the overall factor of X13X24 mentioned above. 

The final step is to impose -D4 x S3 symmetry by summing the quantity (3.3) over all 
dihedral transformations of the external xi,X2,x^,X4 and all permutations of the internal 
xa,xb,xc- Applying this procedure to (3.3), and keeping in mind our choice to normalize 
each distinct term with a factor of 1, leads immediately to the expression given in (2.23). 
It is important to note that because we always impose the -D4 x sl symmetry at the level 
of the integrand, the choice of explicit loop momenta labeling 1 xa, xb, ■ - ■' ha these figures is 
largely irrelevant: they only serve as a representative with which we may present a particular 
numerator explicitly; any other choice would have been acceptable, as it would also generate 
the same rational function after symmetrization. (The labels of loop momenta with only four 
propagators are particularly irrelevant, as they do not appear in the numerator at all, and so 
are never needed to specify an explicit, representative integrand). 

3.2 Dual Conformally Invariant Diagrams 

It is known that there are 1,1,2,10,59 distinct pseudo-conformal diagrams respectively at 
loop-orders L = 1, 2, 3, 4, 5, respectively. However, only 1, 1, 2, 8, 34 of those diagrams actually 
enter the L-loop four-particle integrand with nonzero coefficient; the remaining 2 four-loop 
and 25 five-loop diagrams do not appear in the amplitude. In [62] it was observed that a precise 
characterization can be given which distinguishes the contributing versus non-contributing 
diagrams: the former consist of those pseudo-conformal diagrams which are rendered infrared 
finite when they are evaluated off-shell, i.e. when we take the external momenta pi to satisfy 
pf 7^ 0. (Note that there are additional types of DCI diagrams one can draw when off-shell 
legs are allowed. These have been classified through four-loops in [63], but these additional 
diagrams play no role in our present analysis.) 

To illustrate this point let us consider the simplest five-loop diagram which is pseudo- 
conformal but not DCI, 



X24X 3C 



(3-3) 



X13X 24 



XWX2AX2CX3AX3BX4BX4CXABXACXBC 




whose associated integrand is 



^13^24 



+ sym. (120 terms total) . 



(3.4) 



xieX2CX3DX4bxabxacxadxaexbdxbexcdxce 
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Let us focus on the region of integration where all four of xb,xq,xd and ig simultaneously 
approach xa- Specifically let us suppose that xab,%AC,%AD,%AE are all order 0(p 2 ) as 
p — >• 0. In this limit there are eight vanishing propagators, giving a factor of p ie in the 
denominator, and the integral behaves as 



where the factor p dp in the numerator comes from switching to radial coordinates for 
xb,xc,xd,xe near the point where all of these equal xa- Evidently this integral has a 
divergence near p = 0. 

In contrast none of the diagrams shown in figure 1, nor any of the 8 (34) diagrams which 
are known to contribute to the 4 (5)-loop integrands, suffer from any divergences of this 
sort. We refer to pseudo-conformal diagrams which possess this property as genuinely dual 
conformally invariant (DCI) diagrams. Let us emphasize that this stands as an empirical 
observation, as there is no rigorous understanding of why this peculiar fact should play any 
important role. In particular let us emphasize that DCI diagrams most certainly do have 
infrared singularities when evaluated in any of the familiar regularization schemes such as 
dimensional regularization or the Higgs regulator of [35]. Nevertheless, appreciating the 
utility of this finiteness criterion through five-loops, we adopt the conjecture that it continues 
to hold to all- loop order. 

3.3 Classification of Four-Point DCI Diagrams Through Seven-Loops 

In this section we explain the steps taken to classify all DCI diagrams through seven-loops. 
The classification takes place in two steps: we first used an off-the-shelf program to gener- 
ate all relevant distinct planar topologies, and then wrote our own code to determine, for 
each topology, all independent sets of numerator factors which render the diagram dual con- 
formally invariant. The vast majority of planar graphs (in particular, any which contain a 
triangle — an internal face with only three edges) cannot be rendered DCI with any numerator 
factor, and can be rejected immediately. On the other hand, beginning at five-loops there are 
examples of planar graphs which are rendered dual conformally invariant by more than one 
inequivalent choice of numerator factor (the three occurrences of this phenomenon at five- 
loops are evident in fig. 7). Counting diagrams with different numerator factors separately 
we find 1,1,2,8,34,256,2329 distinct four-point DCI diagrams at one- through seven-loops 
respectively. These come from 1, 1, 2, 8, 30, 197, 1489 distinct underlying topologies (ignoring 
numerators). These results are summarized in table 1. 

The first step in the classification requires the choice of a graph generating program. 
The popular choice qgraf [64, 65], is not well-suited to the present application because of 
our interest exclusively in plane graphs. Note that at this point we must begin to be precise 
about distinguishing a 'planar graph', which is one that admits an embedding onto the plane 
with no self- intersection, from a 1 plane graph', which is a planar graph together with a specific 
choice of plane-embedding. This distinction is important because a given planar graph can 




(3.5) 



- 13 - 



have more than one inequivalent embedding. Two plane graphs which correspond to the 
same underlying planar graph would be considered identical in scalar field theory but must 
be distinguished in gauge theory because they have different color index flow. 

This rather subtle distinction first becomes necessary at six-loops, where we find the 
candidate planar graph shown in fig. 2 which admits two inequivalent plane-embeddings, 
only one of which can be made pseudo-conformal with a suitable numerator factor. Note that 
in this case the pseudo-conformal graph on the left is not DCI (it fails the finiteness criterion 
explained in section 3.2). We strongly suspect that it is not possible for any genuine DCI 
diagram to admit an inequivalent embedding in the plane, but this notion plays no role in 
our analysis and we only mention this phenomenon in passing as a curiosity. 

Returning to our main point: in order to both avoid wasting effort generating non-planar 
graphs in which we are not interested, and to be sure of generating all possible inequivalent 
plane graphs, we turn to the specialized program plantri [66, 67] which is capable of exactly 
what we require. The program has no concept of external lines, so to represent a four-point 
diagram in plantri we first 'fold it up' to a plane graph with (L + 4) faces by tying the 
four external lines together at a new 'vertex at infinity' A given {L + 4)-face plane graph, 
together with the choice of a quadrivalent vertex in the graph, may then be 'unfolded' to a 
planar four-point diagram. The process of folding and unfolding is not bijective since on the 
one hand two inequivalent four-point graphs may become identical when folded up, while on 
the other hand a given plane graph with k > 1 quadrivalent vertices may unfold to fewer than 
k topologically distinct four-point diagrams. The latter happens when a plane graph has a 
discrete symmetry which exchanges two (or more) of its quadrivalent vertices. 

To generate all topologies of interest we 

1. used plantri to enumerate all plane graphs G\,Gz, . . . with vertex connectivity > 1 
and at least one quadrivalent vertex; 

2. unfolded each graph Gi at each of its quadrivalent vertices to obtain various four-point 




Figure 2. Two different plane embeddings of the same planar graph (they are related to each other 
by rotating the bottom edge of the square by 180 degrees out of the plane of the paper while leaving 
the other three edges fixed). The graph on the left can be made pseudo-conformal with a suitable 
numerator factor while the graph on the right cannot. 
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diagrams G\ L> ,G\ Z >,... (taking care to remove any duplicates, as mentioned above); 

3. from the collection of diagrams so generated, we excluded any which have internal 
triangles, since these cannot possibly be made DCI. 

It turns out that the plantri output is most naturally sorted not by L but by the number 
of vertices in the graph. An L-loop DCI graph can have up to 2L + 3 vertices (including the 
vertex at infinity), so in order to classify all diagrams through seven- loops we need all graphs 
with up to 17 vertices. The reader may be interested to know that we found exactly 

1, 0, 1, 1, 6, 13, 71, 308, 1637, 8421, 45229, 243220, 1326433 (3.6) 

'DCI candidate' plane graphs with 5 through 17 vertices respectively. By a DCI candidate 
we mean a plane graph which has at least one quadrivalent vertex with the property that 
unfolding the graph at that point leads to a triangle-free four-point diagram for the denomi- 
nator. 

The four-point diagrams obtained from these DCI candidates were then fed into a custom 
Java program which scanned through all possible inequivalent numerator factors, selecting 
those which make the diagram both pseudo-conformal according to section 3.1 and finite 
according to the criterion reviewed in section 3.2. Many of the DCI candidates do not admit 
any valid numerator, while a few admit more than one inequivalent choice of numerator. 
The final output from our code program was a list of the 1, 1, 2, 8, 34, 256, 2329 distinct DCI 
integrands through seven- loops advertised above; refer also to table 1. 

(These numbers include 1 6-loop and 2 7-loop four-point diagrams with vertex connec- 
tivity 1, which have vertex connectivity 2 when folded up. These diagrams do not contribute 
to scattering amplitudes. The 6-loop example was discussed in figure 12 of [3].) 

As part of this work we in fact generated all DCI candidate plane graphs through 18 
vertices, but found it computationally prohibitive to tackle 19, which would have been nec- 
essary for a complete classification at eight-loops. While an L-loop DCI graph can have as 
many as 2L + 3 vertices, it can on the other hand have as few as L + 4 (examples saturating 
this lower bound exist for L > 3). Therefore, although our classification of DCI diagrams is 
complete only through seven- loops, the data we have amassed contains an enormous number 
of additional DCI diagrams through fourteen- loops, which we have not yet analyzed. 

4 Bootstrapping the Four-Point Integrand Through Seven-Loops 

The way in which the four-point multi-loop integrand can be obtained via the soft-collinear 
bootstrap was illustrated for the three-loop integrand in section 2.4. The condition that 
(logilLi)^ has at most an 0(l/e)-pole in the soft-collinear limit (2.13) is equivalent to the 
condition that the numerator of (log M&p L ' (when all contributing terms are assembled over a 
common denominator) must vanish in this limit (as the denominator is certainly proportional 
to e 2 ). With this in mind, the procedure for obtaining the L-loop integrand M4 from the 
soft-collinear bootstrap can be carried out as follows: 
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# of denominator 


# of distinct 


# of integrands with coefficient: 


Loop 


topologies 


DCI integrands 


+1 


-1 


+2 





1 


1 


1 


1 











2 


1 


1 


1 











3 


2 


2 


2 











4 


8 


8 


6 


2 








5 


30 


34 


23 


11 








6 


197 


256 


129 


99 


1 


27 


7 


1489 


2329 


962 


904 


7 


456 



Table 1. Statistics of four-point DCI integrands through seven-loops. 



1. expand the log of the L-loop amplitude — as illustrated in eqn. (2.5) — and compute 
the contribution of each term involving strictly lower-loop amplitudes to the overall 
numerator in the soft-collinear limit (2.13); call this 'log rest '; 

2. find all L-loop DCI integrands l\ (with each one fully symmetrized with respect to 
permutations of the loop momentum variables A%, . . . , Al, and the L>4 dihedral sym- 
metry of the external momenta), and compute the contribution of each to the overall 
numerator in the soft-collinear limit (2.13); call these contributions (L ) so ft! 

3. log rest and the (L^) S oft are of course still nontrivial polynomials in various x a fe's, and 
our goal is now to determine a collection of numerical constants Cj such that 

iogrcst + ^2ci{l\ L) ) so i t = 0. 

% 

By repeatedly evaluating this equation at sufficiently many random independent values 
of the remaining variables {xa% > • ■ • > X A L as wen as x±, . . . ,X4) this equation can be 
turned into a linear system on the Cj with constant coefficients. The existence and 
uniqueness of the solution, which we find through seven-loops, are the central miracles 
of this paper. 

Our conjecture is that this general strategy extends to all loop-orders. This procedure is 
simple enough to automate, and can be implemented efficiently enough to find all integrands 
through six-loops in less than 2 minutes on a sufficiently powerful computer (given knowledge 
of the L-loop DCI integrands L^). 

For the sake of providing a uniform reference for the reader, the results for four and 
five-loops are reviewed in figures 6, and 7, respectively. For six- and seven-loops, there are 
too many contributions to be practically included graphically here, but they may be viewed 
at [54]. 

Beyond five-loops, there are some interesting features worth noting, which are summa- 
rized in table 1. For example, although all 34 five- loop DCI integrands contribute to the 
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amplitude, this is not generally the case: at six-loops, 27 of the 256 DCI integrands have 
vanishing coefficients; and at seven-loops, 456 of the 2329 DCI have vanishing coefficients. 
The 27 six-loop DCI integrands with vanishing coefficient are illustrated in figure 8. 

More curiously, at six-loops there is the first appearance of a coefficient found to be '+2'; 
and at seven-loops, seven integrands have coefficient +2. These are illustrated in figures 3 
and 4, respectively. It would be fascinating to see if this pattern of coefficients could be 
understood along the lines of the method proposed in [68] . 

One available and relatively simple partial check on our results involves two-particle 
cuts, which have been shown to iterate to all orders in SYM theory [69]. This means that if 
there is an L-loop contribution to the integrand with the property that it can be reduced to a 
product of an L — l-loop contribution and the four-particle tree amplitude by cutting only two 
internal propagators, then that L-loop diagram must inherit the same coefficient as its L — 1- 
loop ancestor. At two- and three-loops all DCI integrands have two-particle cuts; the first 
examples without two-particle cuts are the second and third graphs in fig. 6 (which therefore 
are the first graphs "allowed" by two-particle cut considerations to have coefficients other 
than +1). We have verified the consistency of our results with this principle; at seven-loops 
we find that 814 of the 2329 coefficients correspond to integrals which admit a two-particle 
cut. 

A Appendix 

In this appendix we collect a few useful graph theoretic results about DCI diagrams. 

Lemma 1. Every four-point DCI diagram has an overall factor 0/X13X24. 

Proof. We will show that the overall factor of X13X24 is required in order to pass the finiteness 
criterion of section 3.2. Consider an L-loop pseudo-conformal diagram with P propagators 
(internal edges). By dimensional analysis there must be a total of P — 2L numerator factors. 
Now consider the behavior of the integrand in the limit when all of the x's (both internal 




Figure 3. The only single six-loop integrand contributing with coefficient +2. 



-17- 



Figure 4. The 7 seven-loop contributions with coefficient +2. 




and external) approach a common point y. If (x — y) 2 = 0{p 2 ) for all x's then the integrand 
scales like 



n 2P 



(A.l) 

where we have included the appropriate factor p iL ~ l dp from the integration measure in radial 
coordinates near p = 0. 

Of course the four x's associated to the external are supposed to be fixed, so the fact that 
we have a divergence in this hypothetical limit is not yet of concern. If we remove a single 
external face x a from consideration (that is, hold x a fixed instead of letting it approach y), 
then the propagators adjacent to x a no longer contribute a factor of p~ 2 each to eqn. (A.l), 
but the numerator factors attached to x a also no longer contribute a factor of p +2 . Since 
there are an equal number of each, there is no change in the conclusion and we still have a 
log-divergent integral. 

If we remove two more external faces from consideration this analysis is unchanged unless 
at some step we encounter a numerator factor x™ b connecting two external edges. In order to 
avoid counting the no- longer-contributing factor of p 2m twice we should multiply eqn. (A.l) 
by p 2m , which renders the integral finite for any m > 0. 

Therefore, by the time we've excluded three external faces from the limit (thereby reach- 
ing an honest limit of the integral) , we must have encountered at least one numerator factor 
between two external faces in order to avoid a divergent integral. Since the diagram must 
converge for every possible limit in order to be DCI, we conclude that for any choice of three 
external faces x a , Xj,, x c , two of them must always share a numerator factor (that is, there must 



\ 



Figure 5. A degenerate graph is one in which two external legs are attached to the same vertex. Such 
graphs can be pseudo-conformal (although this example is not), but they cannot be DCI. 



be a positive power of at least one of x a b, Xb c or x ca ). Since x a b = if a and b are adjacent, 
the only way this is possible is for the diagram to contain an overall factor of £13X24- D 

Since each of the four external faces xi, ... ,£4 appears at least once in the numerator 
of every DCI integrand, and conformal invariance requires that each one appears exactly as 
often in the denominator as in the numerator, we have the immediate 

Corollary 1. Degenerate diagrams (see figure 5) cannot be DCI. 

Lemma 2. Every individual four-point DCI diagram diverges as C(l/e 2 ) in the limit (2.13). 

Proof. This is tantamount to proving that every DCI integrand contains a term with both 
xia and X2A m the denominator. A general graph has a denominator of the form 

' - (A.2) 



where each * denotes an arbitrary loop integration variable. Since each external face must 
appear at least once in the denominator, each of the four products in parentheses is non-empty. 

In summing something of the form (A.2) over its Sl x D4 images, the only way to avoid 
ever encountering a factor X1AX2A hi the denominator is for all of the *'s connected to any of 
the four external x a to all be the same, as in the expression 



1 



(xia • • • x\a){x2A • • • x 2 a)(x3A ■ • ■ x 3 a)(x 4 a • • • x a a)x** ■ ■ ■ X** 
However since double poles are not allowed we can have at most 

■^13*^24*^** * * " 



(A.3) 



XIAX2AX3AX4AX** • • • x** 
where the numerator factor is forced on us by Lemma 1 and dual conformal invariance. 



(A.4) 
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Figure 6. The integrands which generate the four-loop amplitude. 




By dimensional analysis the number of x**'s in the numerator must be 2(L — 1) less than 
the number of S in the denominator. Therefore in the limit when all other integration 
variables approach xa (say), the integrand behaves like 

X13X24 d 4 arAP 4(L-1)-1 dp x 13 x 2 a d 4 x A dp 
~ 777 — r\ = (A. 5) 

X 1A X 2 AX3AX4A p^ L - l > X 1A X2AX3 A X4 A p 

which diverges. Therefore we conclude that such a term is diagram is divergent if it does 
not contain a term with both x\ A and x<ia in the denominator. Now, after summing over 
Sl x D4 images, a given diagram will in general contain many terms which each have both x\ A 
and X2A in the denominator. Since there are no possible sources of minus signs these terms 
can never sum to zero identically. This concludes the proof that a DCI diagram is divergent 
unless it contains a term with both x± A and X2A in the denominator. □ 
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Figure 8. The 27 six-loop contributions with vanishing coefficients. Notice that each of the last three 
graphs admit more than one numerator; the coefficient of each possible integrand is separately zero. 
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